Satellite imagery as a management tool for monitoring water clarity across freshwater ponds on Cape Cod, Massachusetts.

Water clarity serves as both an indicator and a regulator of biological function in aquatic systems. Large-scale, consistent water clarity monitoring is needed for informed decision-making. Inland freshwater ponds and lakes across Cape Cod, a 100-km peninsula in Massachusetts, are of particular interest for water clarity monitoring. Secchi disk depth (SDD), a common measure of water clarity, has been measured intermittently for over 200 Cape Cod ponds since 2001. Field-measured SDD data were used to estimate SDD from satellite data, leveraging the NASA/USGS Landsat Program and Copernicus Sentinel-2 mission, spanning 1984 to 2022. Random forest machine learning models were generated to estimate SDD from satellite reflectance data and maximum pond depth. Spearman rank correlations (rs) were "strong" for Landsat 5 and 7 (rs = 0.78 and 0.79), and "very strong" for Landsat 8, 9, and Sentinel-2 (rs = 0.83, 0.86, and 0.80). Mean absolute error also indicated strong predictive capacity, ranging from 0.65 to 1.05 m, while average bias ranged from -0.20 to 0.06 m. Long- and recent short-term changes in satellite-estimated SDD were assessed for 193 ponds, selected based on surface area and the availability of maximum pond depth data. Long-term changes between 1984 and 2022 established a retrospective baseline using the Mann-Kendall test for trend and Theil-Sen slope. Generally, long-term water clarity improved across the Cape; 149 ponds indicated increasing water clarity, and 8 indicated deteriorating water clarity. Recent short-term changes between 2021 and 2022 identified ponds that may benefit from targeted management efforts using the Mann-Whitney U test. Between 2021 and 2022, 96 ponds indicated deteriorations in water clarity, and no ponds improved in water clarity. While the 193 ponds analyzed here constitute only one quarter of Cape Cod ponds, they represent 85% of its freshwater surface area, providing the most spatially and temporally comprehensive assessment of Cape Cod ponds to date. Efforts are focused on Cape Cod, but can be applied to other areas given the availability of local field data. This study defines a framework for monitoring and assessing change in satellite-estimated SDD, which is important for both local and regional management and resource prioritization.

Providing a framework for seagrass mapping in United States coastal ecosystems using high spatial resolution satellite imagery.

Seagrasses have been widely recognized for their ecosystem services, but traditional seagrass monitoring approaches emphasizing ground and aerial observations are costly, time-consuming, and lack standardization across datasets. This study leveraged satellite imagery from Maxar's WorldView-2 and WorldView-3 high spatial resolution, commercial satellite platforms to provide a consistent classification approach for monitoring seagrass at eleven study areas across the continental United States, representing geographically, ecologically, and climatically diverse regions. A single satellite image was selected at each of the eleven study areas to correspond temporally to reference data representing seagrass coverage and was classified into four general classes: land, seagrass, no seagrass, and no data. Satellite-derived seagrass coverage was then compared to reference data using either balanced agreement, the Mann-Whitney U test, or the Kruskal-Wallis test, depending on the format of the reference data used for comparison. Balanced agreement ranged from 58% to 86%, with better agreement between reference- and satellite-indicated seagrass absence (specificity ranged from 88% to 100%) than between reference- and satellite-indicated seagrass presence (sensitivity ranged from 17% to 73%). Results of the Mann-Whitney U and Kruskal-Wallis tests demonstrated that satellite-indicated seagrass percentage cover had moderate to large correlations with reference-indicated seagrass percentage cover, indicative of moderate to strong agreement between datasets. Satellite classification performed best in areas of dense, continuous seagrass compared to areas of sparse, discontinuous seagrass and provided a suitable spatial representation of seagrass distribution within each study area. This study demonstrates that the same methods can be applied across scenes spanning varying seagrass bioregions, atmospheric conditions, and optical water types, which is a significant step toward developing a consistent, operational approach for mapping seagrass coverage at the national and global scales. Accompanying this manuscript are instructional videos describing the processing workflow, including data acquisition, data processing, and satellite image classification. These instructional videos may serve as a management tool to complement field- and aerial-based mapping efforts for monitoring seagrass ecosystems.

Simulated response of St. Joseph Bay, Florida, seagrass meadows and their belowground carbon to anthropogenic and climate impacts.

Seagrass meadows are degraded globally and continue to decline in areal extent due to human pressures and climate change. This study used the bio-optical model GrassLight to explore the impact of climate change and anthropogenic stressors on seagrass extent, leaf area index (LAI) and belowground organic carbon (BGC) in St. Joseph Bay, Florida, using water quality data and remotely-sensed sea surface temperature (SST) from 2002 to 2020. Model predictions were compared with satellite-derived measurements of seagrass extent and shoot density from the Landsat images for the same period. The GrassLight-derived area of potential seagrass habitat ranged from 36.2 km2 to 39.2 km2, averaging 38.0 ± 0.8 km2 compared to an observed seagrass extent of 23.0 ± 3.0 km2 derived from Landsat (range = 17.9-27.4 km2). GrassLight predicted a mean seagrass LAI of 2.7 m2 leaf m-2 seabed, compared to a mean LAI of 1.9 m2 m-2 estimated from Landsat, indicating that seagrass density in St. Joseph Bay may have been below its light-limited ecological potential. Climate and anthropogenic change simulations using GrassLight predicted the impact of changes in temperature, pH, chlorophyll a, chromophoric dissolved organic matter and turbidity on seagrass meadows. Simulations predicted a 2-8% decline in seagrass extent with rising temperatures that was offset by a 3-11% expansion in seagrass extent in response to ocean acidification when compared to present conditions. Simulations of water quality impacts showed that a doubling of turbidity would reduce seagrass extent by 18% and total leaf area by 21%. Combining climate and water quality scenarios showed that ocean acidification may increase seagrass productivity to offset the negative effects of both thermal stress and declining water quality on the seagrasses growing in St. Joseph Bay. This research highlights the importance of considering multiple limiting factors in understanding the effects of environmental change on seagrass ecosystems.

Vertical artifacts in high-resolution WorldView-2 and WorldView-3 satellite imagery of aquatic systems.

Satellite image artefacts are features that appear in an image but not in the original imaged object and can negatively impact the interpretation of satellite data. Vertical artefacts are linear features oriented in the along-track direction of an image system and can present as either banding or striping; banding are features with a consistent width, and striping are features with inconsistent widths. This study used high-resolution data from DigitalGlobe's (now Maxar) WorldView-3 satellite collected at Lake Okeechobee, Florida (FL), on 30 August 2017. This study investigated the impact of vertical artefacts on both at-sensor radiance and a spectral index for an aquatic target as WorldView-3 was primarily designed as a land sensor. At-sensor radiance measured by six of WorldView-3's eight spectral bands exhibited banding, more specifically referred to as non-uniformity, at a width corresponding to the multispectral detector sub-arrays that comprise the WorldView-3 focal plane. At-sensor radiance measured by the remaining two spectral bands, red and near-infrared (NIR) #1, exhibited striping. Striping in these spectral bands can be attributed to their time delay integration (TDI) settings at the time of image acquisition, which were optimized for land. The impact of vertical striping on a spectral index leveraging the red, red edge, and NIR spectral bands-referred to here as the NIR maximum chlorophyll index (MCINIR)-was investigated. Temporally similar imagery from the European Space Agency's Sentinel-3 and Sentinel-2 satellites were used as baseline references of expected chlorophyll values across Lake Okeechobee as neither Sentinel-3 nor Sentinel-2 imagery showed striping. Striping was highly prominent in the MCINIR product generated using WorldView-3 imagery, as noise in the at-sensor radiance exceeded any signal of chlorophyll in the image. Adjusting the image acquisition parameters for future tasking of WorldView-3 or the functionally similar WorldView-2 satellite may alleviate these artefacts. To test this, an additional WorldView-3 image was acquired at Lake Okeechobee, FL, on 26 May 2021 in which the TDI settings and scan line rate were adjusted to improve the signal-to-noise ratio. While some evidence of non-uniformity remained, striping was no longer noticeable in the MCINIR product. Future image tasking over aquatic targets should employ these updated image acquisition parameters. Since the red and NIR #1 spectral bands are critical for inland and coastal water applications, archived images not collected using these updated settings may be limited in their potential for analysis of aquatic variables that require these two spectral bands to derive.

Temporal Stability of Seagrass Extent, Leaf Area, and Carbon Storage in St. Joseph Bay, Florida: a Semi-automated Remote Sensing Analysis.

Seagrasses are globally recognized for their contribution to blue carbon sequestration. However, accurate quantification of their carbon storage capacity remains uncertain due, in part, to an incomplete inventory of global seagrass extent and assessment of its temporal variability. Furthermore, seagrasses are undergoing significant decline globally, which highlights the urgent need to develop change detection techniques applicable to both the scale of loss and the spatial complexity of coastal environments. This study applied a deep learning algorithmto a 30-year time series of Landsat 5 through 8 imagery to quantify seagrass extent, leaf area index (LAI), and belowground organic carbon (BGC) in St. Joseph Bay, Florida, between 1990 and 2020. Consistent with previous field-based observations regarding stability of seagrass extent throughout St. Joseph Bay, there was no temporal trend in seagrass extent (23 ± 3 km2, τ = 0.09, p = 0.59, n = 31), LAI (1.6 ± 0.2, τ = -0.13, p = 0.42, n = 31), or BGC (165 ± 19 g C m-2, τ = - 0.01, p = 0.1, n = 31) over the 30-year study period. There were, however, six brief declines in seagrass extent between the years 2004 and 2019 following tropical cyclones, from which seagrasses recovered rapidly. Fine-scale interannual variability in seagrass extent, LAI, and BGC was unrelated to sea surface temperature or to climate variability associated with the El Niño-Southern Oscillation or the North Atlantic Oscillation. Although our temporal assessment showed that seagrass and its belowground carbon were stable in St. Joseph Bay from 1990 to 2020, forecasts suggest that environmental and climate pressures are ongoing, which highlights the importance of the method and time series presented here as a valuable tool to quantify decadal-scale variability in seagrass dynamics. Perhaps more importantly, our results can serve as a baseline against which we can monitor future change in seagrass communities and their blue carbon.

Acute health effects associated with satellite-determined cyanobacterial blooms in a drinking water source in Massachusetts.

BACKGROUND: The occurrence of cyanobacterial blooms in freshwater presents a threat to human health. However, epidemiological studies on the association between cyanobacterial blooms in drinking water sources and human health outcomes are scarce. The objective of this study was to evaluate if cyanobacterial blooms were associated with increased emergency room visits for gastrointestinal (GI), respiratory and dermal illnesses. METHODS: Satellite-derived cyanobacteria cell concentrations were estimated in the source of drinking water for the Greater Boston area, during 2008-2011. Daily counts of hospital emergency room visits for GI, respiratory and dermal illnesses among drinking water recipients were obtained from an administrative record database. A two-stage model was used to analyze time-series data for an association between cyanobacterial blooms and the occurrence of illnesses. At the first stage, predictive autoregressive generalized additive models for Poisson-distributed outcomes were fitted to daily illness count data and daily predictive variables. At the second stage, residuals from the first stage models were regressed against lagged categorized cyanobacteria concentration estimates. RESULTS: The highest cyanobacteria concentration (above the 75th percentile) was associated with an additional 4.3 cases of respiratory illness (95% confidence interval: 0.7, 8.0, p = 0.02, n = 268) compared to cyanobacteria concentrations below the 50th percentile in a two-day lag. There were no significant associations between satellite derived cyanobacterial concentrations and lagged data on GI or dermal illnesses. CONCLUSION: The study demonstrated a significant positive association between satellite-derived cyanobacteria concentrations in source water and respiratory illness occurring 2 days later. Future studies will require direct measures of cyanotoxins and health effects associated with exposure to cyanobacteria-impacted drinking water sources.

Assessing cyanobacterial frequency and abundance at surface waters near drinking water intakes across the United States.

This study presents the first large-scale assessment of cyanobacterial frequency and abundance of surface water near drinking water intakes across the United States. Public water systems serve drinking water to nearly 90% of the United States population. Cyanobacteria and their toxins may degrade the quality of finished drinking water and can lead to negative health consequences. Satellite imagery can serve as a cost-effective and consistent monitoring technique for surface cyanobacterial blooms in source waters and can provide drinking water treatment operators information for managing their systems. This study uses satellite imagery from the European Space Agency's Ocean and Land Colour Instrument (OLCI) spanning June 2016 through April 2020. At 300-m spatial resolution, OLCI imagery can be used to monitor cyanobacteria in 685 drinking water sources across 285 lakes in 44 states, referred to here as resolvable drinking water sources. First, a subset of satellite data was compared to a subset of responses (n = 84) submitted as part of the U.S. Environmental Protection Agency's fourth Unregulated Contaminant Monitoring Rule (UCMR 4). These UCMR 4 qualitative responses included visual observations of algal bloom presence and absence near drinking water intakes from March 2018 through November 2019. Overall agreement between satellite imagery and UCMR 4 qualitative responses was 94% with a Kappa coefficient of 0.70. Next, temporal frequency of cyanobacterial blooms at all resolvable drinking water sources was assessed. In 2019, bloom frequency averaged 2% and peaked at 100%, where 100% indicated a bloom was always present at the source waters when satellite imagery was available. Monthly cyanobacterial abundances were used to assess short-term trends across all resolvable drinking water sources and effect size was computed to provide insight on the number of years of data that must be obtained to increase confidence in an observed change. Generally, 2016 through 2020 was an insufficient time period for confidently observing changes at these source waters; on average, a decade of satellite imagery would be required for observed environmental trends to outweigh variability in the data. However, five source waters did demonstrate a sustained short-term trend, with one increasing in cyanobacterial abundance from June 2016 to April 2020 and four decreasing.

Satellite remote sensing to assess cyanobacterial bloom frequency across the United States at multiple spatial scales.

Cyanobacterial blooms can have negative effects on human health and local ecosystems. Field monitoring of cyanobacterial blooms can be costly, but satellite remote sensing has shown utility for more efficient spatial and temporal monitoring across the United States. Here, satellite imagery was used to assess the annual frequency of surface cyanobacterial blooms, defined for each satellite pixel as the percentage of images for that pixel throughout the year exhibiting detectable cyanobacteria. Cyanobacterial frequency was assessed across 2,196 large lakes in 46 states across the continental United States (CONUS) using imagery from the European Space Agency's Ocean and Land Colour Instrument for the years 2017 through 2019. In 2019, across all satellite pixels considered, annual bloom frequency had a median value of 4% and a maximum value of 100%, the latter indicating that for those satellite pixels, a cyanobacterial bloom was detected by the satellite sensor for every satellite image considered. In addition to annual pixel-scale cyanobacterial frequency, results were summarized at the lake- and state-scales by averaging annual pixel-scale results across each lake and state. For 2019, average annual lake-scale frequencies also had a maximum value of 100%, and Oregon and Ohio had the highest average annual state-scale frequencies at 65% and 52%. Pixel-scale frequency results can assist in identifying portions of a lake that are more prone to cyanobacterial blooms, while lake- and state-scale frequency results can assist in the prioritization of sampling resources and mitigation efforts. Satellite imagery is limited by the presence of snow and ice, as imagery collected in these conditions are quality flagged and discarded. Thus, annual bloom frequencies within nine climate regions were investigated to determine whether missing data biased results in climate regions more prone to snow and ice, given that their annual summaries would be weighted toward the summer months when cyanobacterial blooms tend to occur. Results were unbiased by the time period selected in most climate regions, but a large bias was observed for the Northwest Rockies and Plains climate region. Moderate biases were observed for the Ohio Valley and the Southeast climate regions. Finally, a clustering analysis was used to identify areas of high and low cyanobacterial frequency across CONUS based on average annual lake-scale cyanobacterial frequencies for 2019. Several clusters were identified that transcended state, watershed, and eco-regional boundaries. Combined with additional data, results from the clustering analysis may offer insight regarding large-scale drivers of cyanobacterial blooms.

Quantifying national and regional cyanobacterial occurrence in US lakes using satellite remote sensing.

Cyanobacterial harmful algal blooms are the most common form of harmful algal blooms in freshwater systems throughout the world. However, in situ sampling of cyanobacteria in inland lakes is limited both spatially and temporally. Satellite data has proven to be an effective tool to monitor cyanobacteria in freshwater lakes across the United States. This study uses data from the European Space Agency Envisat MEdium Resolution Imaging Spectrometer and the Sentinel-3 Ocean and Land Color Instrument to provide a national overview of the percentage of lakes experiencing a cyanobacterial bloom on a weekly basis for 2008-2011, 2017, and 2018. A total of 2321 lakes across the contiguous United States were included in the analysis. We examined four different thresholds to define when a waterbody is classified as experiencing a bloom. Across these four thresholds, we explored variability in bloom percentage with changes in seasonality and lake size. As a validation of algorithm performance, we analyzed the agreement between satellite observations and previously established ecological patterns, although data availability in the wintertime limited these comparisons on a year-round basis. Changes in cyanobacterial bloom percentage at the national scale followed the well-known temporal pattern of freshwater blooms. The percentage of lakes experiencing a bloom increased throughout the year, reached a maximum in fall, and decreased through the winter. Wintertime data, particularly in northern regions, were consistently limited due to snow and ice cover. With the exception of the Southeast and South, regional patterns mimicked patterns found at the national scale. The Southeast and South exhibited an unexpected pattern as cyanobacterial bloom percentage reached a maximum in the winter rather than the summer. Lake Jesup in Florida was used as a case study to validate this observed pattern against field observations of chlorophyll a. Results from this research establish a baseline of annual occurrence of cyanobacterial blooms in inland lakes across the United States. In addition, methods presented in this study can be tailored to fit the specific requirements of an individual system or region.

Performance across WorldView-2 and RapidEye for reproducible seagrass mapping.

Satellite remote sensing offers an effective remedy to challenges in ground-based and aerial mapping that have previously impeded quantitative assessments of global seagrass extent. Commercial satellite platforms offer fine spatial resolution, an important consideration in patchy seagrass ecosystems. Currently, no consistent protocol exists for image processing of commercial data, limiting reproducibility and comparison across space and time. Additionally, the radiometric performance of commercial satellite sensors has not been assessed against the dark and variable targets characteristic of coastal waters. This study compared data products derived from two commercial satellites: DigitalGlobe's WorldView-2 and Planet's RapidEye. A single scene from each platform was obtained at St. Joseph Bay in Florida, USA, corresponding to a November 2010 field campaign. A reproducible processing regime was developed to transform imagery from basic products, as delivered from each company, into analysis-ready data usable for various scientific applications. Satellite-derived surface reflectances were compared against field measurements. WorldView-2 imagery exhibited high disagreement in the coastal blue and blue spectral bands, chronically overpredicting. RapidEye exhibited better agreement than WorldView-2, but overpredicted slightly across all spectral bands. A deep convolutional neural network was used to classify imagery into deep water, land, submerged sand, seagrass, and intertidal classes. Classification results were compared to seagrass maps derived from photointerpreted aerial imagery. This study offers the first radiometric assessment of WorldView-2 and RapidEye over a coastal system, revealing inherent calibration issues in shorter wavelengths of WorldView-2. Both platforms demonstrated as much as 97% agreement with aerial estimates, despite differing resolutions. Thus, calibration issues in WorldView-2 did not appear to interfere with classification accuracy, but could be problematic if estimating biomass. The image processing routine developed here offers a reproducible workflow for WorldView-2 and RapidEye imagery, which was tested in two additional coastal systems. This approach may become platform independent as more sensors become available.